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Abstract 

On the basis of the concept of the growing role of nonadiabatic effects of the non- 
conservation of the quantum number K, a theory has been developed of the phe- 
nomenon which has been given the name of backbending. Above the transition point, 
for J > Jc, all the values — J < K < J are equally probable. An investigation is made 
of the singularities possessed by the ordering parameter (proportional to the spectro- 
scopic quadrupole moment of a nonspherical nucleus), the rotational angular velocity 
and the moment of inertia of a nucleus at the Curie point. Formulas have been derived 
for the intensity of quadrupole radiation in the more symmetric n-phase J > Jc- By 
analyzing the experimental values of the moments of inertia belonging to the n-phase, 
the radius of the mass distribution in the nucleus was determined. It agrees with the 
radius of the proton distribution derived from data on the scattering of electrons by 
nuclei. On the basis of the simplest form of the singularity of the parametric derivative 
of the Hamiltonian of the system a general theory of zero-temperature second-order 
phase transitions is developed in the Appendix. 

Introduction. Estimate of the critical value Jc of the 
rotational quantum number. Rotational density ma- 
trix 



In recent decades there perhaps has not been in nuclear physics an event more outstanding 
than the discovery made at the beginning of the 1950's by A. Bohr and B. Mottelson of 
rotational levels in nonspherical nuclei. Having confirmed as a subsidiary result the well- 
known assertion of quantum mechanics concerning the inability of a perfectly spherical body 
to undergo purely mechanical rotation (cf., for example, Refs. [1,2]), this discovery became 
the starting point for a detailed study of the so-called rotational bands in the energy spectra 
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of nonspherical nuclei. The initial simplest theoretical treatment is expressed by the widely 
known formula 

E{J)^^AJ+i) (1) 

for the energy levels. Here J is the rotational quantum number, i.e., the total angular 
momentum of the whole system (nuclear spin); /' is the "adiabatic" moment of inertia which 
refers to the neighborhood of the origin of the band and the value of which is determined 
from experiment. Avoiding as far as possible complications which do not have particular 
significance as to principle, we shall in future as a rule have in mind the case of a rotational 
band based on the ground state of an even-even nucleus. Then in fact the allowable values 
are J = 0,2,4, . . . 

Progress of the experiment upward along the band stimulated an ever more critical at- 
titude to formula (1). One of the methods of generalizing it consists of the following: we 
adopt the point of view which appears to be natural that the conserved quantity J is the 
only significant physical (thermodynamic) characteristic of which the energy is a function 
within the band under consideration in the case of not too light a nucleus. Any scalar con- 
structed from the components of the vector J will be reduced to the combination J{J + 1). 
Therefore 

oo 

£^m(J)= J]C,[J(J+1)]^ (2) 

s=l 

(the reason for indexing the rotational energy of the nucleus by m will become clear from 
subsequent discussion). Prom a practical point of view doubt immediately arises whether 
neglecting terms of the series (2) after a finite their number would lead to a situation that 
is sufficiently stable from a purely computational point of view. However with time an even 
deeper difficulty involving matters of principle became apparent. About six years ago due 
to the application of new experimental methodology penetration began into the domain of 
still greater values of J [3-6] . This led to the discovery of a very characteristic phenomenon 
which was called backbending ^. When a noticeable effect was observed it consisted of the 
following: the principal characteristics 

dE (fE dVL ff 

(Q, is the rotational angular velocity, / is the moment of inertia) of the rotational motion 
of a nucleus undergo a sharp change over a very narrow region of the rotational band. The 
angular velocity of rotation Vl here falls off precipitously in spite of an increase in the total 
angular momentum J. It is also very characteristic that in some cases immediately beyond 
the transition point the rotational spectrum turns out to be close to an equidistant one. Por 
example, quite typical in this regard is the rotational band based on the ground state of the 
nuchde W^^*^. The segment J = 10 — 14 in this case contains only two intervals. But over 
this segment the angular velocity Ml (measured in energy units) manages to decrease by 
~ 100 keV — by approximately 30%. The sector J = 12 — 16 which is minimally displaced 

-'^This is due to the external appearance of curves that illustrate the position of rotational levels of the 
band and constructed in terms of somewhat artificially chosen coordinates. In the present work we shall not 
make use of such coordinates. 
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in the positive direction also contains two energy intervals. But the second of them exceeds 
in the first one by only 3%, and the increase in angular velocity of rotation here is just as 
insignificant (cf., formula (3)). 

At first glance it might appear that the discovery of the phenomenon of backbending 
compromises the very basis for arguing in favor of relations of the type (1) or (2). However 
it appears to us to be not quite logical and excessively hasty to cast doubt on the determining 
role of the independent (and, with macroscopic accuracy, essentially unique) thermodynam- 
ical variable J. In particular, the arguments leading to formula (2) by no means exclude 
the possibility that this series would have a finite radius of convergence. In accordance with 
the spirit of experimental data the point J — Jc oi the intersection of the boundary of the 
circle of convergence with the real axis could then have an important physical meaning. If 
nuclear states situated on opposite sides of this point differ by some qualitative property 
of a symmetry type, then for J > Jc the energy E{J) must be a different function which 
does not coincide with the analytic continuation of expression (2). Due to the presence of a 
singularity at J = Jc such an analytic continuation could, strictly speaking, turn out to be 
not single-valued and even complex. 

In order to discuss the specific physical nature of the phenomenon we turn to the prop- 
erties of the component J ■ n = (n is a unit vector along the axis of the non-spherical 
nucleus) of the total angular momentum of the nucleus along its symmetry axis. Although 
this is not a matter of principle and does not in any manner affect exact calculations, in order 
to make the subsequent discussion more easily visualizablc one can, if one so wishes, call 
upon a crudely classical vector model; one can set = {Jz)max = J and then K = Jcos6. 
The quantum number K is not conserved since the operator J ■ n does not commute with 
the Hamiltonian of the system. However for J <^ near the origin of the rotational band 
the commutator is small. It is well known that this corresponds to adiabatic slowness of 
rotation. As a result we obtain 

K^Ko = const, wk = Skko, J < Jc, (4) 

where wk is the probability of having the corresponding value of K in the quantum state 
under consideration of the nucleus as a whole. Here formula (1) is valid. In practice we most 
frequently have to deal with the simplest case Kq = for even-even nuclei. 

But as one moves along the rotational band the admixture of components with K ^ Kq 
increases more and more rapidly. A rough estimate making use of perturbation theory leads 
to a relation of the type 

^/;k(x [2|X|/(J- |i^|)]^-l^le^-l^l 

to express the law of growth of a given component near the threshold J = \K\ where it first 
appears. High above the threshold we obtain the expression Wk{J) oc J^'^', which remains 
sufficiently steep and also favors the appearance of large K J. As long as perturbation 
theory is applicable the initiating component K — Kq appears to serve, in agreement with 
unitarity, as a source of growth for the others. However, when the latter cease being relatively 
small, then all the component become qualitatively on the same footing in the process of 
their, so to say, interaction with each other. The result will be a kind of equilibrium, an 
equidistribution with respect to K {wk = 1/(2 J -|- 1) = const), which will be attained, say, 
at the point J — Jc- Of course, here one should not expect any appreciably sharp change 
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in the state of the system as such, for example a discontinuity in its energy. Physically the 
important feature of such a new situation is contained in a different attribute — in the change 
in the symmetry of the rotational state of the nucleus. It can be easily understood that an 
equidistribution with respect to K corresponds to isotropy, to an equal probability of all 
spatial directions of the vector n (we recall that we are dealing with an individual specific 
quantum state of the nucleus as a whole: J and are fixed quantities!). At the same 
time in the Hamiltonian of the system as a whole one also should not be able to discern 
any favored directions for the vector n in empty space. It would be naive to think that 
this correspondence will be violated in the course of a further increase in J when physically 
there are even fewer grounds for some values of the quantum number K to become favored 
compared to others. Having acquired stability the more symmetric isotropic phase will not 
lose it for J > Jc- 

The attainment of isotropy with respect to the directions of n can also be interpreted from 
a somewhat different point of view. In a sufficiently strong rotational "field" the mechanical 
angular momenta of the individual quasiparticles aline parallel to the vector J and cease 
to aline along the vector n. Therefore the latter turns out to be "free," i.e., it is in fact 
distributed isotropically (for J > Jc). This constitutes the difi^erence from the old adiabatic 
region J <^ where the internal state of the nucleus is determined, roughly speaking, by 
its deformation and, in turn, determines the approximately conserved value = Jn ~ Xq. 
At the transition point itself the number of aligned quasiparticle is of order unity, i.e., 

Je ~ i ~ fc/i? > 1 (5) 

{kf is the limiting momentum of the Fermi distribution; R is the nuclear radius). This 
estimate does not contradict experimental data (cf., also Sec. 5). 

Thus, J = Jc is the point of a phase transition of the second kind — the zero-temperature 
Curie point. In analogy with the terminology adopted in Ref. 7 we shall refer to the states 
of a nonsphcrical nucleus with a completely isotropic distribution of its axis in space as 
the n-phase, and we shall call the region J < Jc oi "spontaneous symmetry violation" the 
m-phase. We then have 

WK = ^—, K^^^— y ■^(■^+^) , j>j^ (6) 

2J+1' 3' 2J+l^j 3 ' - ^ ^ 

{9 and (p are the spherical angles defining the direction of the vector n with respect to 
fixed axes). We emphasize that these properties of the n-phase presuppose a most essential 
indeterminancy in the value of the polar angle 6 which does not in any way diminish with 
increasing J. In the m-phase the relations (6) are violated as a result of the "ordering" of 
the rotational state. Having in mind, primarily, the special case Kq = and taking into 
account the fact that in the adiabatic region J <^ Jc (cf., formulas (1) and (4)) a change in 
the quantum number J has little effect on the quantity K'^, we shall, in order to be specific, 
set 

X;^^X^<^^^^, J<Jc. (7) 

k=-J 

It can be easily verified that for a state which is a superposition of components with 
different values of K it is not possible to construct any purely rotational wave function 
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which depends only on n. For J ^ Jc the quantum mechanical description of the rotation 
of a nonspherical nucleus is achieved with the aid of the density matrix p(n, n'). However, 
due to the conservation of the 2;-component of the total angular momentum the "azimuthal" 
wave function 

0.,* = (27r)-VV^^ M = J,. (8) 

nevertheless exists. In the quasiclassical case = {Jz)max = J ^ 1 it corresponds to the 
so-called regular precession of an ordinary symmetric top (cf., for example, Ref. 8). As 
regards the polar angle 9^ sufficiently far above the Curie point the density matrix which 
depends upon it must have the following properties: p{9^ 9') = for 9 ^ 9', while the angular 
dependence of the diagonal element of p{9, 9') corresponds to isotropy in three-dimensional 
space. However the formula appears more attractive if one utilizes a representation which is 
based not on angles, but on the more customary quantum numbers J, M and K. Then we 
have ^ 

Pj^jSJ, M, K- J', M\ K') = 5jj^5j,j,5mj.5m'jMk' (9) 

in a state of angular momentum Jq and its component J^. For diagonal elements the region 
of applicability of this relation is wider and encompasses the whole n-phasc including the 
Curie point. But for J + 0, apparently, non-diagonal (with respect to K) elements of 

the density matrix which reflect the correlations between different K also become signiflcant. 
It is even more difficult to form judgements concerning the specific form of the rotational 
density matrix in the m- phase J < Jc- For J = Jc it has some kind of a singularity which 
hinders the formal continuation of the corresponding expression into the "foreign" region 
J > Jc- The physically important functional of the rotational density matrix — the so-called 
order parameter — has a radical singularity as J Jc — (cf., next section). 



2 The average ("spectroscopic") quadrupole moment 
as an order parameter. Discontinuity in the rota- 
tional angular velocity of the nucleus. 

If we keep (6) and (7) in mind, the following definition suggests itself 

3^2 



J(J + 1) 



^l-3cos2^ (10) 



for the order parameter which characterizes deviation of the m-phase from isotropy (the 
approximate expression involving the angle refers to the quasiclassical case Jz — J ^ 1 
and is stated here only for ease of visualization). In the important special case Kq = it 
in fact varies from unity to zero at the Curie point. Naturally, for a given J and rj quite 
different nuclear states are in principle conceivable — both rotational states, and ones differing 
microscopically. However, we shall assume that for J < Jc one can utilize the concept of the 
energy of that one of these states which is "most nearly in equilibrium" and which we shall 
denote by E{J,ri). 

Intending to analyze the situation for sufficiently small Jc — J in the neighborhood of 
the Curie point we expand the general expression for the energy in powers of 7] and limit 
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ourselves to three terms: 



E{J, ri) = En{J) - A{J)ri + lD{J)r]^ (11) 

The function En{J) in fact refers to the unordered n-phase. The inequahties 

A(J) > 0, D(J) ^ D(J^) = D>0 (12) 

express the advantages of positive rj and the stabihty of the energy minimum with respect to 
this quantity. Substitution into (11) of the true value of rj determined from the equilibrium 
condition 

transforms E{J,r]) into the energy E^{J) of the rotational levels of the m-phase 

V=^, En{J)-E^iJ) = \^. (14) 

Figure 1, which shows the curve for the energies of the phases, requires some explanation. 
At the "end-point" J = Jc the function Em{J) has a certain singularity. However, the 
question is: Which of the derivatives of the function will be first affected by this, i.e. which 
of them will become infinite at this point? According to (3) this cannot happen in the case of 
the first derivative: there are no slightest theoretical or experimental bases for supposing that 
the rotational velocity of nuclei becomes infinite. Moreover, as may be seen from formula (3), 
the second derivative also does not become infinite. A zero value of the moment of inertia 
agrees neither with the fact that the nucleus is non-spherical, nor with the experimentally 
observed tendency for moments of inertia in the m-phase to vary (cf. Sec. 5). Thus, although 
the dotted line extrapolation of the Ejn{J) curve is somewhat arbitrary, at least the first two 
derivatives have meaning along it. As a result of this assumption illustrated by Fig. 1 that 
at J = Jc a simple point of intersection 

En - EmO^ Jc - J (15) 

of the curves under investigation occurs is not contradictory and agrees well with the nature 
of the experimental data. 

Comparison of (15) with (14) yields 

A{J) = a{Jc - J)^/^ a > 0. (16) 

Consequently 

r]^a{Jc- Jf'/D (17) 
and ^ 

K-£^m = ^^(Jc- J). (18) 

Differentiating relation (18) we obtain a discontinuity in the rotational angular velocity of 
the nucleus at the point of "backbending" 

A(;ifi) = a^/D. (19) 
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Here A{hfl) = hflmc ~ ^nc- As far as one can judge according to experimental data, the 
thermodynamic inequahty 

is never violated. 

It is well known [1,2] that for the value averaged over the state = J of the 2;-component 
of the quadrupole moment (the so-called spectroscopic quadrupole moment) of a nonspherical 
nucleus the following formula is valid 

(J+l)(2J + 3) ^^^^ 

Here, Qo is the component of the "collective" macroscopic quadrupole moment along the nu- 
clear axis wholly determined by its axially-symmetric deformation. Taking (10) into account 
we have 

'? = -2JT3'=^°''- '''' 

Now comparing this with (17) we verify that both the quadrupole moment Q and also the 
order parameter vanish at the Curie point according to (Jc — J)^/^. In the n-phase J > Jc 
there can no longer be any collective (spectroscopic) quadrupole moment because of its 
isotropy. 



3 Position of the rotational levels in the n-phase J > Jc. 
Moment of inertia. 

First of all we touch in a couple of words upon the problem of the stabihty of either of the 
phases. The corresponding condition 

has a fairly natural appearance ^. The symmetric n-phase becomes unrealizable for J < 
because of the absolute instability of the state in which the mechanical angular momenta of 
the quasiparticles aligning along the vector J ignore the direction of the vector n determined 
by the deformation (cf., the in a certain sense opposite nature of the states which are in fact 
realized in the adiabatic region J <C Jci cf., also the Introduction). Therefore we have 



= (24) 

J=Jc 

^It is not difRcult to verify that negative moments of inertia correspond to absolute instability both 
thermodynamic and purely mechanical. The latter in the present case has the following meaning; free 
rotation of a body corresponds, as is well known, to thermodynamic equilibrium and is not accompanied by 
friction (cf., for example, Refs. 9 and 10). If the regular precession with respect to the angle cj) is treated in 
accordance with the laws of classical mechanics (cf. the Introduction, formula (8) and the text referring to 
it), then for negative moments of inertia the principle of least action is violated. 
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and 

'^^J-J. (25) 
near the Curie point. Substitution into the second of formulas (3) yields 

I^-r—r, J-Jc<^Y- (26) 

J — Jc Jo 

Here j is a certain constant coefficient; 

lo = ^Mi?2 (27) 

is the solid body value of the moment of inertia; M = rUnA is the nuclear mass; m„ is the 
nucleon mass; 

R = roAV3 (28) 

is the nuclear radius. We arrive at an important conclusion that the moment of inertia of a 
nonspherical nucleus becomes infinite in accordance with the law (26) if one approaches the 
Curie point from above. This is what explains the observed approximate uniform spacing 
of levels in corresponding segments of the rotational spectra (cf., the Introduction). In 
the opposite limiting case of large J — Jc, for a sufficiently clearly pronounced lining-up of 
mechanical angular momenta of the quasiparticles along the direction of the vector J the 
moment of inertia of the Fermi system becomes the same as that for a solid body: 

/^/o, J-Jc>f. (29) 

-'0 

If the right-hand sides of (26) and (29) are added, then an interpolation formula will be 
obtained for the moment of inertia of the n-phase, which satisfies both limiting cases. By 
itself it, of course, is devoid of any to some extent deep physical foundations. Moreover, in 
the case of a literal acceptance of this interpolation a situation would arise of a monotonic 
approach to the solid body limit Iq from above, while experimentally, apparently, a minimum 
was observed in the variation of the moment of inertia (cf.. Sec. 5). However, in a purely 
technical respect the interpolation formula is a convenient route for obtaining correct limiting 
expressions for fln{J) a-nd En{J) by means of integrating it twice. Retaining everywhere in 
addition to the terms containing constants of integration, a single term of the corresponding 
expansion we obtain the following formulas 

m^mnc + —{J - Jcf, 



for J — Jc <^ j / Iq and 



6j 



^2 

Ml = hflnc + —{J - Jc), 
-'0 

E^Eo + Klnc{J -Jc) + ^{J- Jc)^ 



(31) 
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for J-J^:^ j/Io. 

It is desirable to keep in mind that, generally speaking, the last formula predicts the 
position of the rotational levels with good relative but not absolute accuracy. We note that 
for nuclides for which j / Iq ~ 1 the domain of applicability of formulas (31) is broadened 
encompassing practically the whole n-phase. 



4 Electric quadrupole radiation in the n-phase 

The basic property of the n-phase can be formulated as the absence of an observable (average) 
quadrupole moment; 

(Qik) = 0. (32) 

This does not mean that its matrix elements non-diagonal with respect to the levels of the 
band also vanish. They are responsible for the quadrupole radiation the probability of which 
we can write in the form 

-"-|^(«^»- (33) 

We have for the sake of brevity symbolically denoted by ((Q^)) (reduced intensity of the 
transition) the square of the absolute value of the matrix element between the final state 
J — 2 and the initial state J summed over the final orientations of the nuclear spin. The 
transition frequency of interest to us amounts to c<j = 2Q. 

For the subsequent discussion it is convenient to express this collective quadrupole mo- 
ment in terms of the components of the vector n: 

Qik = ^Qoi^riiUk - Sik). (34) 

Prom this it can be seen that it is immaterial over which state the square of the tensor is 
averaged: 

Ql = IQI- (35) 
With the aid of the usual rules for matrix multiplication we verify the validity of the relation 

T.\iQ^>^)a\" = IQI (36) 
6 

Here the indices a and b enumerate the individual states of a nonspherical nucleus irrespec- 
tively of to which rotational band they happen to belong. 

Above the Curie point the quadrupole radiation can be calculated quasiclassically. Setting 
(j) = fit we calculate in accordance with (34) the classically varying quadrupole moment. 
Then in accordance with the isotropy of the n-phase it is averaged uniformly over cos^ 
{K — J cos 9 in the quasiclassical case). Finally the square of the absolute value of the 
spectral component a; = 2Q is substituted in formula (33) in place of ((Q^)). As a result of 
elementary calculations we obtain ^ 

{{Q')) = \Ql ^=^^^0, J>Jc. (37) 



Since the operator for collective quadrupole moment is diagonal with respect to K, the correlations 
mentioned at end of the Introduction do not affect the result. 
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Table 1 



Nucleus 




MeV 


MeV 


Nucleus 


Jc 


ibi ij fnc 

MeV 


A{h^lmc) 

MeV 




13 


530 


200 




15 


345 


70 




11 


410 


175 




17 


370 


120 




13 


490 


235 




>15 


310 


>20 




11 


600 


365 


Ybl64 


15 


350 


110 




>17 


340 


>25 




15 


340 


90 




13 


275 


75 


Ybl68 


15-17 


330 


~30 




>15 


310 


>15 


Y^iyo 


17 


365 


60 




>15 


320 


>45 




15 


340 


110 


J3yi56 


>17 


370 


>40 




15-17 


300 


~5 




13 


275 


95 


tttI 79 


>17 


360 


>40 


I3yl58 


15 


335 


45 




13 


320 


105 


Qyl60 


15 


310 


~25 




>15 


350 


>40 




13 


380 


120 




15 


340 


100 




15 


360 


125 




>13 


385 


>30 




15 


345 


70 




>15 


460 


>210 



Thus, within the bounds of the n-phase all the radiative transitions have the same reduced 
intensity. In addition to the transition responsible for the real radiation there also exists 
an analogous transition in the opposite direction upwards along the band. The sum of the 
reduced intensities of both transitions is obtained by doubling the right hand side of the 
first of formulas (37). Comparing this with the sum rule (36) we find that it is saturated 
to one-third within the band under consideration, while the remaining two-thirds refer to 
transitions to other rotational bands. 

5 Comparison with experiment 

Data on the states of nonsphcrical even-even nuclei with high spins were taken primarily 
from Ref. 11. The variation of the rotational angular velocity hQm in the m-phase was ex- 
trapolated graphically into a relatively small region near the Curie point (cf., also comments 
on Fig. 1 in Sec. 2), while the latter was identified with the position of the energy interval 
which already belongs, as far as the available data enable us to judge, to the n-phase. The 
results of this treatment of data arc shown in Table 1. The two cases marked by an asterisk 
refer to the rotational band originating from the 0^ level of the excitation of f3 oscillations 
(the fully symmetric oscillation of a nonsphcrical nucleus). The quite unsystematic varia- 
tions in the magnitude of the increment A{hQ) of the rotational velocity from one nucleus 
to the next are noteworthy. The critical angular velocity hft^c behaves in a much more 
stable manner, but it also is subject to random fluctuations. It is somewhat more difficult 
to judge how real on the average is the tendency to a certain amount of decrease in Hflmc 
with increasing atomic weight. 

In addition to those given in Table 1 there are six more nuclides in which rotational levels 
have been discovered up to J^ax = 14 - 18: W^^^ Wl^^ Os^^^, Os^^^ Os^^^, and Th"^^^. No 
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Table 2 





Cei30 




J3yl58 


Eri62 


Ybl64 


Ybi66 




08^82 


Jc 

J/Io 


11.0 
2.11 


12.4 
2.12 


14.8 
6.76 


15.7 
2.63 


13.0 
3.98 


15.5 
2.11 


13.6 
1.88 


11.9 
4.20 



significant effect was observed. However, let us turn to the rotational velocities attained. 
For the three osmium isotopes mentioned above we have HQ = 280 — 315 keV. This is lower 
than the values characteristic of the end of Table 1 and is close to the lower limit of critical 
velocities in the whole tabic. For both tungsten isotopes hQmax = 300 keV. The thorium 
nucleus in the experiment attained the value hQ = 200 keV. A conclusion suggests itself that 
in the nuchdes "suspected" of not showing the effect there was simply a lack of sufficient 
rotational velocity to equalize the intensities of the components of the wave function with 
different values of K. Therefore there is every reason to suppose that with a further increase 
in J the situation will become clarified and, as usual, a phase transition will be found. 
We note, in order to avoid misunderstanding, that even in the case of a vanishingly small 
increment A{hn) (a situation close to this already occurred in the case of Hf^™, cf.. Table 1) 
it still will be possible to recognize the X-phase by its moment of inertia. 

The study of the variation of the moment of inertia of a nonspherical nucleus in both 
phases is of considerable interest. The most characteristic and complete data in this respect 
are those concerning the principal rotational band for W^^''; cf.. Fig. 2. The increase in the 
moment of inertia in the m-phase can be qualitatively understood in the spirit of the Le 
Chatelier-Brown principle as a "resistance to an external action". Indeed, the rotational 
perturbation responsible for the growth in the components with K ^ Kq is determined not 
directly by the mechanical angular momentum J, but by the angular velocity Qm- Therefore 
the system tries to diminish the derivative d{hfl)/dJ. ^ On the other side of the Curie point 
in the m-phase, the moment of inertia drops somewhat below Iq, passes through a minimum 
and then rapidly tends to approach the solid body asymptote from below. Analogous non- 
monotonic behavior in the variation of the quantity (PEn/dJ"^ was also observed for the 
nuclide Yb^®^. An interesting case was encountered in the /9-vibrational band of dysprosium 
^j^yieo* according to the notation of Table 1). From all indications here the very first two 
experimental points belonging to the n-phase turned out to be situated on opposite sides of 
the minimum in the moment of inertia. A final verification of the hypothesis will have to be 
postponed until the position of the rotational levels of the band with J > 18 is established. 

On the whole the data concerning the m-phase are so far relatively scarce. When both 
experimental values closest to the Curie point belong to the left wing of the moment of 
inertia curve formula (26) yields 

j J2-J1 ' " (hyih-ihyi), ■ ^^^^ 

These relations enable us to obtain the "extrapolated" (from the domain of the n-phase) 

*An analogous interpretation can also be given to the increase in the moment of inertia of the n-phase on 
the right wing of its curve (cf., Fig. 2) — only here the system resists the alignment which begins to acquire 
a macroscopic scale of the mechanical angular momenta of the individual quasiparticles along the direction 
of the vector J (cf., the Introduction). This is what explains the existence of the minimum in the variation 
of the moment of inertia of the n-phase which will be discussed subsequently. 
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position of the point of phase transition, and also to evaluate j / Iq- For the eight nuclides 
the results are shown in Table 2. Basically good agreement is observed with values of Jc 
shown in Table 1 determined by a different method. The greatest discrepancy occurs in the 
case of osmium. However, it must be stated that the value of Jc = 12 given in Table 2 falls 
just in the middle of that segment of the rotational band where a decrease in the rotational 
angular velocity occurs. The region occupied by the phase transition is here anomalously 
wide and encompasses three energy intervals. Possibly this is in some way associated with 
the closeness of osmium to the point of phase transition of nonspherical nuclei into spherical 
ones [7] (cf., also the striking anomalies at the beginning and at the end of Table 1). We 
also note that for specific nuclides the values of the ratio j / Iq vary no less randomly than 
the increments A(/ifi). 

The universal character of the asymptotic properties of the n-phase far from the Curie 
point in principle opens a path to the determination of the radius of the mass distribution 
in a nucleus. Unfortunately at the present time such high spins of rotational states have not 
yet been attained for which the relation (29) could be confidently regarded as a sufficiently 
strict equality. However one can try to circumvent this difficulty in the following manner. 
On the right wing of the moment of inertia curve for the n-phase, where it is increasing, the 
corresponding values are lower than the rigid body values. This very fact provides us with 
a lower bound on the possible values of the nuclear radius. A less rigorous but still quite 
plausible assumption consists of the following: in cases of the type selected for Table 2 the 
experimental values of the moment of inertia of the n-phase closest to the Curie point are 
higher than the rigid body values. This provides an upper bound on the radius. Making 
use of the fact that data are available concerning eight different nuclides we narrow down as 
much as possible the bounds of the inequality: 

1.08F < ro < l.llF (39) 

(1F= 10~^^ cm). Thus, with reasonable accuracy we have 

ro = I.IF. (40) 

This agrees with the measurements of the radius for the distribution of electrical charge in 
the nucleus from electron scattering [12]. In the preceding discussion the value (40) was 
utilized everywhere in the calculations of the rigid body moment of inertia using formulas 
(27) and (28). 

For a comparison of formulas (30) and (31) with experiment one must have in the n-phase 
a sufficiently large number of experimental points referring to the same nuclide. The most 
favorable case is the one of W^^°. Since the value of j/ Iq given in Table 2 is not large we may 
use formula (31). It is not difficult to understand the reason for the good agreement over 
the whole n-phase demonstrated in Table 3. Being quite well founded under the condition 
J — Jc S> j / Iq formula (31) at the same time has a correct interpolation behavior: as J Jc 
it gives an experimental value fine for the rotational angular velocity of the nucleus which 
is known to be correct (cf., also the remark at the end of Sec. 3). In order to illustrate the 
application of formula (30) we in contrast select a nuclide with the greatest value of j/Io', in 
accordance with Table 2 this will be Dy^^^. As Table 3 shows, the theory agrees well with 
experiment. 
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Table 3 



W^^" Nucleus 


Dy^^s Nucleus 


J 


Ml (keV) experim. 


hfl (keV) theory 


J 


Ml (keV) experim. 


hfl (keV) theory 


13 


216.4 




15 


289.0 




15 


223.2 


240 


17 


295.5 


295.7 


17 


265.0 


273 


19 


312.9 


313.3 


19 


308.3 


306 


21 


339.0 


342.0 


21 


342.6 


339 









Table 4 



Nucleus 


Jc 


h^lmc (keV) 


A{hn) (kcv) 


Ho^^^ 


>37/2 


355 


>100 




>37/2 


340 


>70 


Ho^ei 


>37/2 


325 


>55 




39/2 


380 


~25 


Ybi65 


35/2 


305 


70 


Lu^e^ 


>27/2 


300 


>15 



We now pass to the question of the electric quadrupole radiation from a non-adiabatically 
rapidly rotating nucleus. In a recently published experimental paper [13] a study was made 
of the radiation within the hmits of the principal rotational band of Ce^^^. Two observed 
transitions belong to the n-phase. For comparison with previous formulas we reproduce the 
value oi ioT K — Kq — 0: 

(cf., for example, Ref. [14]; wc neglect the corrections ~ 1/ J^). Dividing the right hand side 
of the first of formulas (37) by (41) we obtain the reduced intensity F measured relative to 
its purely adiabatic value: 



Setting J = 12 we have Ftheor — 0.48. According to Ref. 13, F^xp — 0.65 ± 0.13 for 
the transition 12"*" — > lO"*", while the next transition 14'^ — > 12"*" has the intensity F^^p — 
0.64 ±0.28. 



In conclusion we touch upon the case of odd nuclei. Because it is difficult to sort out the 
numerous rotational bands until now there are relatively few data concerning them. We know 
that the principal source of experimental information concerning the problem of interest to 
us are the quadrupole 7 quanta corresponding to transitions with AJ = 2. Therefore the 
usual adiabatic "rotational band" characterized by certain Kq and parity should be separated 
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into two appropriate bands and treated as distinct bands. ^ Results referring to six nuclides 
and based on the experimental papers [15-17] are given in Table 4. They give no cause for 
doubting the universal character of the phenomenon. 

We express our gratitude to A.I. Baz', I.I. Gurevich, I.M. Pavhchenkov and K.A. Ter- 
Martirosyan for discussions of questions touched upon in this paper and its results. The 
authors are particularly grateful to I.M. Pavhchenkov for many stimulating discussions of 
the nature of experimental data concerning the phenomenon of backbending. 



Appendix 

It appears to be of some interest to analyze certain specific features of zero-temperature 
phase transitions of the second kind from a more general point of view. 

Let the energy E of the ground state of a macroscopic body and the state itself 
depend on a certain thermodynamic quantity x, and the corresponding condition for stability 
be represented for the sake of brevity and simplicity in the form 

analogous to (23). We shall assume that for x < Xc a. spontaneous breaking of symmetry 
occurs which is characterized by the ordering parameter rj. An analytic continuation of the 
symmetric n-phase into the region x < Xc is not physically realizable because of 



d'E 



dx"^ 



0. (A.2) 



Formally this can be only the result of some kind of singularity of the paramedic deriva- 
tives of the Hamiltonian H{x) acting on the variables ^ of the system. In accordance with 

H{x)%{C;x)^Ei{x)%iC;x) (A.3) 

it determines the eigenfunctions and the energy eigenvalues (the Latin subscript enumerates 
the stationary states of the system which parametrically depend on x). The simplest imagin- 
able possibility is the following: the singularity in the derivative dH/ dx leads to the "break" 
in the curve of E{x) of the type shown in Fig. 1. Then proceeding in the same manner as in 
Sec. 2, i.e., actually replacing J in formulas (11) and (14)-(18) by x we obtain the following 
results: 

T) = a{xc - x)/D, 

En — Em — ^^(-^c — x)^ (A. 4) 



dEjn 



dx 



dE„ 



dx 



^Thc division into two "sub-bands" is equally well-founded also from a purely theoretical point of view. 
Indeed, one does not expect any considerable differences in principle from even-even nuclei in the present 
case, while it is desirable to exclude the not entirely clear and non-macroscopic effect of the "alignment" of 
the odd nucleon. It is well known that for ii'o = 5 a similar effect occurs even in the adiabatic region where 
it is more easily susceptible of being investigated (cf., for example, Ref. [1]). 
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However, we emphasize that so simple a possibihty exists only when the following important 
condition is satisfied: the thermodynamic quantity dE/dx\,^^,^^ must not reduce to a func- 
tional of the wave function ^'((^;a;o). In the opposite case a discontinuity in the functional 
can only be the result of a discontinuous variation of the function itself. But this would 
already be a first-order transition and not a second-order phase transition in which we are 
interested for which we always have 

*m(C;3^c) = *n(e;2;c)- (A.5) 

We consider the case when the quantity x has the properties of a dynamic variable (this is 
possible — cf., for example, Ref. 7 — although not necessarily so). In the presence of a velocity 
X an additional kinetic energy proportional to its square appears, and the total energy can 
be expressed in the form 

8 = lB{x)x^ + E{x). (A.6) 

Here B{x) is a mass coefficient. Suppose that the system is moving in the negative direction 
along the x axis. At first sight the hypothesis seems plausible that it, generally speaking, 
will approach the point Xc with a certain finite velocity But then in the relation 

r] = jx (A. 7) 

the condition 

7 = 0, x>x^ (A.8) 

must hold also for x = Xc- the domain of the more symmetric n-phase extends up to the point 
of phase transition inclusively [9] . Having all this in mind we apply the law of conservation 
of energy directly at the Curie point: 

In accordance with the second of formulas (A. 4) (cf., also (A.5)), we ignore the "potential 
energy" E{x). The system cannot pass to the left with a finite velocity x^ — differentiation of 
the first of formulas (A.4) shows that in violation of condition (A.8) Bit x — Xc the coefficient 
7(3;) would change discontinuously from zero to infinity. Thus, i:^ = within the framework 
of the foregoing assumption and, consequently, 

Bnixc) = 0. (A.IO) 

Switching on adiabatically smoothly at i — > 00 the variation x{t) of the coordinate and 
solving the time-dependent Schrodinger equation with the Hamiltonian H{x{t)) we obtain 
the well-known expression 



Bn{x) = 2h''J2 



Ei — E, 



n 



dx 



(A.ll) 



for the mass coefficient (in order to ascribe to the notation a more concrete character we 
denote by the index n the ground states i — within the domain of the n-phase x > Xc)- 
Since all the denominators are positive substitution into (A.IO) yields 



^)=0, x^x.. (A.12) 
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Regarding the amplitude it as a result of normalization turns out to be purely 

imaginary and can be everywhere made equal to zero with an appropriate choice of the 
indefinite phase multiplier ex]){ia{x) in the wave function of the ground state. Then all the 
components of the derivative d^n/ dx along the basis vectors of the Hilbert space turn out 
to have zero values at the Curie point 

d^n/dx\^^^^ = 0. (A.13) 

Expanding the parametric derivative d'^n/dx in series and now measuring the coordinate x 
from the position of the point of phase transition we represent the mass coefficient of the 
n-phase in the form 

Bn{x) = 86x^ (A.14) 



where 



Ei — En 



dx 



(A.16) 



4 ^ - 

For the solution of equations of classical mechanics it is simplest of all to make use of their 
first integral (A. 6): 

m- / i£/bY^\±ty^' a x>o 

"'^^^ - \ ±{2£/B„,y/H if x<Q. ^^-^^^ 

Here the time is so reckoned that at the instant t = the system would be passing through 
the point a; = of the phase transition and the energy S is measured from the "potential 
energy" E{0). The upper sign corresponds to motion in the positive direction of the x 
axis and the lower sign to motion in the negative direction. By Bm we have denoted the 
limiting value of the mass coefficient of the m-phase as x — > —0. The singularity which the 
solution of (A.16) in fact has as re +0 deprives the proportionality relationship (A. 7) of 
a strict mathematical interpretation directly at the origin. By this the paradox is removed 
concerning the simultaneous retention both of the law of conservation of energy and the 
nature of the phase transition as the system passes through the Curie point. 

Strictly speaking the states that have been considered arc not entirely equilibrium states, 
since the macroscopic motion is accompanied by inevitable friction. In the general case in 
order to achieve complete equilibrium it is necessary to add to the Hamiltonian a term of the 
form —Xx which does not affect the nature of the phase transition. Equilibrium at x — xq 
corresponds to A = dE / dx\^^^^. If the zero-temperature second-order phase transitions have 
to a large extent a common nature, then the result (A.13) must not depend on whether the 
quantity x docs in fact have the properties of an autonomous dynamical variable. Indeed, 
it is sufficient to require complete thermodynamic equilibrium — "rest", in order that the 
abstract capability in principle of the system to be "in motion" would in no way manifest 
itself physically. Prom this point of view the condition (A.13) imposed on the wave function 
of the n-phasc at the Curie point appears to be of a sufficiently general nature. 

We now return to the dynamic case. The solution (A.16) was classical. The question 
arises whether the limiting expansion (A.14) for the variation of the mass coefficient of the 
n-phase might not lose its concrete physical content when quantum effects are taken into 
account. In order to check this, we turn to the uncertainty relation 

A£Ax - hx (A.17) 
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(cf., for example, Ref. [1]). According to (A. 16) we have 

X = ±{s/by/y2x, 




(A.18) 



From this it can be seen that the requirement Ax <^ x cannot be satisfied only for 



Thus, there exists a possibility in principle to check by a purely classical method the validity 
of the expression (A. 14) down to the smallest values of x. 

In the foregoing friction accompanying motion was not taken into account explicitly. 
Generally speaking, the question concerning dissipative processes is quite complicated and 
requires a more concrete investigation in each individual case. However, in the region where 
the n-phase is sufficiently close to the point of phase transition the situation is simplified. 
Indeed, projecting the time-dependent Schrodingcr equation onto the base vector system 
defined by (A. 3), we establish that the connection between the amplitudes of states with 
different i is realized exclusively by means of coefficients which in addition to the velocity 
X contain matrix elements of the form {'ik\d'ii/dx). They must also be responsible for the 
irreversible quantum jumps from level to level which are the cause of dissipation. But, as 
is clear from (A. 13), such matrix elements vanish at the Curie point and near it friction 
is suppressed. Therefore results which concern the limiting behavior of the n-phase in the 
neighborhood of the phase transition point remain valid. 

Finally, we briefiy touch upon the relation of the foregoing to the special case x = J 
of nuclear rotation considered in the present paper. The complete wave function of a non- 
spherical nucleus can be written as a product 



By ^ we have denoted variables commuting with J and M and with one another (among 
their number is contained also the macroscopic quantity K = J ■ n). The wave functions 
^j(C) depend on J as on a parameter. They obey an equation of the form (A. 3) in which 
the role of the "Hamiltonian" is played by the block diagonal with respect to J and M of 
the complete Hamiltonian for the nucleus (parametrically dependent on J), while rotational 
energy levels serve as eigenvalues. Since the quantity HQq = dE/dJ\j^j^^ docs not reduce to 
a functional of ^'jolO the rotational velocity fl can undergo a discontinuity at the transition 
point as the energy itself E{J) undergoes a continuous variation along the band, and this 
is in fact observed. The discontinuous diminution of Q signifies in essence a change in the 
"coupling scheme" of angular momenta in the system (cf.. Introduction). 
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